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Scientific  Progress 


1 .  Objectives  and  approach 
This  research  is  two-fold  with  the  following  focuses; 

(1)  Development  of  systematic  approaches  for  calibration  of  subsurface  transport  models  using  hard  and  soft  data  on  system 
characteristics  and  behavior  (the  original  research  focus,  Tasks  1-5  below) 

(2)  Temporal  and  spatial  distribution  of  soil  moisture  in  the  heterogeneous  vadose  zone  as  affected  by  near-surface 
atmospheric  boundary  conditions  (additional  work  related  to  landmine  detection,  Tasks  6-8  below) 

This  first  focus  extends  and  combines  recently-developed  tools  for  parameter  estimation  and  geostatistical  simulation  for  use 
with  groundwater  flow  and  transport  models  used  by  the  Army  in  field  investigations  involving  groundwater  quality  and  quantity. 
The  second  task  focuses  on  improving  our  understanding  and  quantification  of  the  distribution  of  soil  moisture  in  the  vadose 
zone  as  affected  by  land-atmospheric  interaction,  evaporation  from  soils,  evapotranspiration  and  land  mine  detection.  The 
primary  aims  of  the  project  were  to: 

•  Develop  an  improved  understanding  of  how  to  add  in  a  systematic  way  the  geostatistical  description  of  aquifer 
heterogeneity  to  the  calibration  process. 

•  Develop  methods  to  condition  parameter  estimation  to  hard  and  soft  data. 

•  Test  our  ability  to  use  sparsely  sampled  data  to  determine  calibration  parameters  at  the  large  scale  where  simulations  are 
conducted. 

•  Extend  the  findings/methods  developed  above  for  more  complex  problems  involving  multi-phase  flow.  More  specifically,  we 
will  investigate  spatial/temporal  soil  moisture  behavior  in  the  shallow  subsurface  (including  around  buried  objects)  under 
various  conditions  such  as  water  table  fluctuation/heat-induced  evaporation  etc.  To  accomplish  this,  a  data  set  has  to  be 
generated  under  well-controlled  conditions.  Focusing  on  understanding  and  quantification  of  the  distribution  of  soil  moisture  in 
the  vadose  zone  (air-water  two-phase  flow  problem),  we  added  new  tasks  (Tasks  6,  7,  and  8). 

•  The  data  sets  generated  in  the  highly  complex  system  but  under  well-controlled  conditions  can  also  be  used  other 
modelers  to  test  their  inversion  schemes. 

Our  approach  uses  existing  forward  and  inverse  modeling  codes  and  new  experimental  data  on  head  response  to  various 
pumping  excitations  and  solute  spills  generated  in  a  three-dimensional  test  bed  facility  at  the  Center  for  Experimental  Study  of 
Subsurface  Environmental  Processes  (CESEP)  at  the  Colorado  School  of  Mines.  The  approach  we  use  includes  the  following 
tasks,  some  of  which  have  been  revised/added: 

Task  1:  Forward  modeling  code  selection 

A  number  of  computer  codes  that  have  the  capability  to  simulate  flow,  reactive  transport  and  mass  transfer  exist.  The  goal 
of  this  research  task  was  to  review  existing  models  that  were  used  by  DOD  and  other  government  agencies  for  decision  making 
and  regulatory  action. 

Task  2:  Review  and  evaluation  of  inverse  codes 

A  primary  tool  that  was  used  in  this  research  was  an  inversion  model  that  uses  data  on  excitation  and  the  corresponding 
responses  to  estimate  the  system  parameters.  The  goal  of  this  research  task  was  to  select  the  appropriate  inversion  code 
through  initial  testing  with  currently  available  data. 

Task  3:  Review  of  methods  used  in  soft  data  collection  and  synthesis 

In  the  original  proposal  this  research  task  was  proposed  to  gather  information  on  the  state  or  practice  of  model  calibration 
at  DOD  sites.  Discussions  with  WES  collaborators,  it  was  determined  that  the  model  calibration  still  relies  on  the  use  of  both 
soft  and  hard  data  and  inversion  is  not  routinely  used.  Hence,  the  efforts  under  this  task  focused  more  on  investigating  the 
effect  of  quantity  of  hard  and  soft  data  on  the  parameter  estimation  accuracy  (see  Task  5). 

Information  such  as  degree  to  which  formal  inverse  modeling  approaches  were  used,  steps  in  the  calibration  process,  utilization 
of  prior  information,  criteria  used  in  determining  good  match,  types  of  soft  data  used  and  degree  of  importance  attached  to  such 
data,  how  issues  of  uncertainty  were  handled  etc.  will  be  gathered  through  a  survey.  In  this  work,  data  was  categorized  as 
either  hard  data  or  soft  data,  depending  on  the  degree  of  quantitative  uncertainty  included  in  the  data. 

Hard  Data:  The  data  collected  at  a  site  on  the  past  behavior  or  from  controlled  large-scale  tests.  The  past  behavior  is  recorded 
through  observations  of  head  or  drawdown  in  wells  and  solute  concentrations  in  monitoring  wells.  Controlled  large-scale  tests 
may  include  pumping  tests  and  natural  gradient  tracer  tests. 

Soft  Data:  All  other  information  that  is  either  qualitative  in  nature  or  are  collected  at  discrete  points  or  localized  areas  in  the 
system.  The  qualitative  data  includes,  geology,  large-scale  stratigraphy  and  geophysical  test  results.  Point  data  include  data 
collected  using  soil  cores  and  cone  penetrometers.  Localized  data  includes  results  from  slug  tests,  push-pull  tests  and  forced 
gradient  tracer  tests. 


Task  4:  Experimental  design 

A  major  component  of  this  research  project  involved  the  design  and  use  of  experimental  data  generated  in  the  laboratory  to 
develop  and  validate  calibration  processes  that  systematically  use  both  hard  and  soft  data.  The  goal  of  this  task  was  to 
generate  a  comprehensive  data  set  under  controlled  conditions  in  intermediate  scale  test  tanks.  The  key  steps  in  the  proposed 
method  were: 

•  Pack  a  three-dimensional  tank  with  a  defined  heterogeneity.  Both  structured  and  stochastically  generated  random  packing 
configurations  are  used.  Well-characterized  test  sands  with  accurately  known  hydraulic  conductivities,  dispersivity,  porosity, 
reactive  properties  and  retention  characteristics  are  used. 

•  Flow  and/or  transport  simulations  are  conducted  to  generate  hard  data  on  system  response  to  various  system  excitations. 
These  simulations  correspond  to  water  table  drawdown  due  to  pumping,  plume  generation  due  to  a  spill,  dissolution  from  a 
NAPL  spill,  etc. 

•  Use  the  generated  data  to  conduct  automatic  calibration  using  the  appropriate  simulator  and  inversion  code.  As  the  system 
characteristics  are  known  exactly  from  the  designed  packing,  it  is  possible  to  determine  the  accuracy  of  the  model  calibration 
and  uncertainty  of  predictions. 

•  Soft  data  were  introduced  in  a  systematic  manner  to  evaluate  how  the  calibration  got  improved  and  uncertainty  of 
predictions  got  reduced. 

Task  5:  Method  and  protocol  development 

This  task  involved  the  analysis  of  results  from  the  previous  tasks  to  develop  a  systematic  procedure  to  calibrate  models  using 
both  hard  and  soft  data.  The  goal  was  to  use  the  knowledge  gained  through  analysis  conducted  using  the  experimental  data  to 
come  out  with  a  calibration  protocol  with  validated  guidelines.  After  reviewing  the  results  from  Tasks  1  through  4,  it  was  inferred 
that  the  effect  of  quantity  of  hard  and  soft  data  on  the  parameter  estimation  accuracy  had  to  be  further  investigated  in  detail, 
and  we  decided  not  to  pursue  the  development  of  a  general  calibration  protocol  (based  only  on  the  data  set  that  we  have 
generated)  that  we  originally  proposed  in  the  proposal.  Therefore,  in  the  new  Task  5,  we  investigated  value  of  soft  and  hard 
data  for  improving  model  calibration  accuracy  (e.g.,  Sakaki  et  al.  [2008],  PDF  attached). 

Task  6:  A  data  set  on  soil  moisture  distribution  in  heterogeneous  vadose  zone 

An  understanding  and  quantification  of  the  distribution  of  soil  moisture  in  the  vadose  zone  is  important  in  many  applications 
involving  land-atmospheric  interaction,  evaporation  from  soils,  evapotranspiration  and  land  mine  detection.  Natural  soil 
heterogeneity  in  combination  with  the  flux  conditions  at  the  soil  surface  creates  complex  spatial  and  temporal  distribution  of  soil 
moisture  in  the  vadose  zone.  Capillary  barrier  effects  at  the  interfaces  of  variation  soil  textures  contribute  to  this  complexity.  We 
added  the  following  tasks  and  generated  data  sets  on  soil  moisture  distribution;  Task  6-1)  in  a  simple  two-dimensional 
heterogeneous  system  under  drainage/wetting  conditions  (e.g.,  Limsuwat  et  al.  [2008],  PDF  attached),  Task  6-2)  in  simple 
two-dimensional  heterogeneous  systems  subjected  to  evaporation  at  surface  (e.g.,  Sakaki  et  al.  [2008],  PDF  attached),  and 
Task  6-3)  in  three-dimensional  heterogeneous  vadose  zone  (e.g.,  Smits  et  al.  [2008],  PDF  attached).  The  goal  of  this  research 
task  was  to  develop  a  data  set  that  describes  the  spatial  and  temporal  distribution  of  soil  moisture  in  simple  heterogeneous 
systems  and  in  three-dimension  a  spatially  correlated  random  field  under  various  boundary  conditions  (e.g,  water  table 
fluctuation,  evaporation  at  the  soil  surface). 

Task  7:  Evaporation  from  ID  column:  Experiment  and  modeling 

The  goal  of  this  task  was  to  perform  controlled  experiments  under  transient  conditions  of  soil  moisture  and  temperature 
using  soil  with  accurately  known  hydraulic/thermal  properties  and  use  this  data  to  test  existing  theories  and  develop  appropriate 
numerical  models.  Water  vapor  flow  under  varying  temperature  gradients  was  implemented  based  on  a  concept  that  allows 
non-equilibrium  liquid/gas  phase  change  with  gas  phase  vapor  diffusion.  In  order  to  validate  this  new  approach,  we  developed  a 
long  column  apparatus  equipped  with  a  network  of  sensors  and  generated  data  under  well-controlled  thermal  boundary 
conditions  at  the  soil  surface.  Water  saturation,  capillary  pressure,  temperature,  relative  humidity  and  column  weight  to  record 
total  mass  of  water  in  the  column  were  continuously  monitored.  Results  from  numerical  simulations  based  on  the  conventional 
equilibrium  and  non-equilibrium  approaches  were  compared  with  experimental  data.  This  knowledge  is  applicable  to  many 
current  environmental  problems  to  include  the  simulation  of  contaminant  transport  and  volatilization  in  the  shallow  subsurface 
and  water  content  fluctuation  in  the  vadose  zone  interacting  with  the  atmosphere.  For  more  details,  see  Smits  et  al. 

Task  8:  2D  tank  experiment  and  modeling  on  moisture/temperature  behavior  around  a  buried  object 

The  objectives  of  this  work  were;  1)  to  develop  a  unique  large  tank  apparatus  equipped  with  a  network  of  recent  sensor 
technologies  for  automated  and  continuous  monitoring  of  moisture,  temperature,  and  relative  humidity,  2)  to  generate  precision 
data  under  well-controlled  transient  heat  boundary  conditions  to  test  a  2D  heat  and  mass  transfer  model  for  soil  containing  a 
buried  object  over  a  diurnal  cycle,  3)  to  compare  model  results  and  experimental  data  and  4)  to  expand  our  exploration  of  the 
effects  of  the  soil  environment  on  landmine  detection  by  investigating  the  influence  of  the  buried  objects  themselves,  soil 


texture,  and  surface  conditions  on  water  content  and  soil  temperatures  above  and  below  various  objects  buried  at  10cm  depth 
and  away  from  it.  A  major  finding  of  this  study  is  that  the  thermal  signature  of  a  buried  object  strongly  depends  on  the  complex 
interaction  between  water  content,  soil  thermal  and  hydraulic  properties  and  geographical  location.  Temperature  differences 
were,  with  some  exceptions,  larger  in  wet  than  in  the  dry  sands,  which  suggests  that  soil  water  may  help  improve  thermal 
signatures  for  some  objects.  Thermal  and  saturation  contrasts  were  generally  very  different  for  a  buried  landmine  than  for  other 
buried  objects.  It  is  our  hope  that  this  data  and  numerical  model  will  improve  the  comprehension  and  interpretation  of  sensor 
imagery  and  lead  to  the  development  of  more  robust  signal  processing  techniques.  More  details  will  be  available  in  Smits  et  al. 
[manuscript  to  be  submitted  soon] 


2.  Transfer  of  data/knowledge 

ARO-USACE  Engineer  Research  and  Development  Center  (ERDC)  supported  researchers  at  Colorado  School  of  Mines  (CSM) 
have  provided  the  ERDC  with  experimental  data  sets  on  the  1)  calibration  of  subsurface  transport  models  using  hard  and  soft 
data  on  system  characteristics  and  behavior  and  2)  temporal  and  spatial  distribution  of  soil  moisture/temperature  in  the  vadose 
zone  that  can  be  used  to  better  understand  fundamental  processes  and  validate  numerical  models. 

The  Pis  have  been  in  contact  with  Dr.  Stacy  Howington  and  Jackie  Hallberg  at  the  WES  to  discuss  how  the  data  and 
research  findings  are  to  be  used  by  the  modeling  group.  A  group  that  includes  Dr.  Stacy  Howington  and  a  research  associate 
Ms.  Jackie  Hallberg  visited  the  CESEP  at  CSM  lab  to  discuss  the  data  that  have  been  collected  and  any  additional  sampling 
that  might  have  to  be  done  to  meet  their  needs.  We  transferred  the  data  we  have  generated  to  them.  Ms.  Hallberg  is  working  on 
parameter  estimation  for  WES’s  countermine  work  and  for  their  other  large  simulators.  In  this  research  she  has  coupled  their 
favorite  flow/transport  model  (ADH)  with  PEST  and  is  ready  for  a  test.  Dr.  Howington  also  indicated  Dr.  Tim  Kelley’s  interest  in 
the  dataset  as  a  test  of  his  optimization  code.  Dr.  Howington  has  provided  additional  funding  for  this  task. 

The  ERDC  has  been  involved  over  the  past  5  years  in  increasing  landmine  detection  and  reducing  false  alarms  through  the 
testing  of  automated  target  recognition  software  (ATRs)  on  synthetic  IR  imagery.  This  work  has  been  funded  through  the 
Countermine,  Near  Surface  Phenomenology,  and  GEOTACS  programs.  The  work  has  focused  on  generating  synthetic  IR 
images  which  are  added  as  test  cases  for  ATR  software.  The  production  of  synthetic  imagery  allows  for  a  variety  of  scenarios 
to  be  tested  that  would  be  costly  or  difficult  to  test  in  the  field.  Because  landmine  sensors  exploit  soil  and  environmental 
conditions  to  discern  between  mines  and  other  objects,  most  mine  detection  technologies  require  that  the  spatial  and  temporal 
variability  of  key  environmental  conditions  such  as  climate,  vegetation,  soil  type,  depth  of  ground  water  table,  and  topography 
be  understood.  If  these  factors  and  the  ability  to  model  them  in  a  variety  of  domains  become  well  defined,  then  sensor  and 
algorithm  simulations  can  more  realistically  be  tailored  to  particular  operational  scenarios  and  technologies.  The  experimental 
data  sets  provided  by  CSM  provide  a  systematic  experimental  study  of  the  effects  of  the  near  surface  boundary  conditions  on 
soil  moisture  and  temperature  distributions  in  the  shallow  subsurface.  Increasing  our  knowledge  of  the  effects  of 
geohydrologic/thermal  properties  and  behaviors  on  the  landmine  signature  will  ultimately  help  to  properly  interpret  sensor 
imagery.  Experimental  results  can  be  used  to  validate  numerical  results  from  the  ERDC’s  Computational  TestBed. 
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